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Abstract 

We study the posibility of large isospin fluctuations in high-energy heavy-ion 
coUisions by assuming that pions are produced semiclassically both directly and in 
pairs through the isovector channel. 

The leading-particle effect and the factorization property of the scattering am- 
plitude in the impact parameter space are used to deflne the classical pion fleld. In 
terms of the joint probability function Pji^{no, nj) for producing no neutral and n_ 
negative pions from a deflnite isospin state //s of the incoming leading-particle sys- 
tem we calculate the two pion correlation parameters „ and the average number 
of neutral pions ((?^o)n_) as a function of negative pions (n_) produced. We show 
that only direct production of pions without isovector pairs leads to large isospin 
fluctuations. 
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1 Introduction 



In recent years several cosmic-ray experiments [1] have reported evidence for the ex- 
istence of Centauro events characterized by an anomalously large number of charged 
pions in comparison with the number of neutral pions, indicating that there should 
exist a strong long-range correlation between two types of the pions. The negative 
results of the accelerator searches for the Centauro events at CERN (2,3,4) have 
suggested that threshold for their production must be larger then 900 GeV. 

Such long-range correlations are possible if pions are produced semiclassicaly 
and constrained by global conservation of isospin [5-10]. 

Although the actual dynamical mechanism of the production of a classical pion 
field in the course of a high-energy collision is not known, there exist a number of 
interesting theoretical speculations [11-16] that localized regions of misahgned chiral 
vacuum might occur in ultrahigh-energy hadronic and heavy-ion collisions. These 
regions become coherent sources of a classical pion field. In early models [5,6], 
however, the coherent production of pions was taken for granted and considered as 
a dominant mechanism. 

These models also predict strong negative correlations between the number 
of neutral and charged pions. In fact, the exact conservation of isospin in a pion 
uncorrelated jet model is known [7,8] to give the same pattern of neutral/charged 
fluctuations as observed in Centauro events. This strong negative neutral-charged 
correlation is believed to be a general property of the direct pion emission in which 
the cluster formation ( or the short-range correlation between pions ) is not taken 
into account [9,10,17,18,19]. 

In this paper, following the approach of our earlier paper (18) we consider the 
leading-particle effect as a source of a classical pion field in the impact parameter 
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space. Pions are assumed to be produced from a definite isospin state of the incoming 
two-particle system both directly and in isovector pairs. 

Coherently emitted isovector clusters decay subsequently into pions outside the 
region of tttt interaction. We discuss the behavior of the joint probability distribution 
Pi 13(^0, n_) of neutral and charged pions, the variation of the two-pion correlation 
functions /^^^ and the average number of neutral pions {{no)n_) as a function of the 
number of negative pions (n_) produced. 

We support the conclusion of (18,19) that the large isospin fluctuations are 
consequence of a direct production of pions. 
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2 The Eikonal S matrix with Isospin 



At high energies most of the pions are produced in the central region. To isolate 
the central production, we adopt high-energy longitudinally dominated kinematics, 
with two leading particles retaining a large fraction of their incident momenta. We 
assume that the collision energy is large enough so that the central region is free of 
baryons. 

The basic assumption of the independent pion-emission model, neglecting the 
isospin for a moment, is the factorization of the scattering amplitude T„(s, 6; 1 . . . n) 
in the b space: 

^ ^ -n-l n 

Tn{s, 6; 1 . . . n) = 2sf{s, b)^l[ J{s, b; q,), (1) 
where, owing to unitarity, 

|/(s,6)P=e-^(^'^) (2) 

and 

n{s,b)^ Jdq\J{s,b;q)\' (3) 

denotes the average number of emitted pions at a given impact parameter b. The 
function | J(s, b; q) p, after the integration over b, controls the shape of the single- 
particle inclusive distribution. A suitable choice of this function also guarantees that 
the energy and the momentum are conserved on the average during the collision. 

The inclusion of isospin in this model is straightforward [6] . The factorization 
of Tn in the form (1) is a consequence of the pion field satisfying the equation of 
motion 

{a + ij^^)n{s,b;x) =j{s,b;x), (4) 
where j is a classical source related to J(s, b; q) via the Fourier transform 

J{s,b;q)^ I d''xe'^^j{s,b;x). (5) 
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The standard solution of Eq.(4) is usually given in terms of in- and out-fields 
that are connected by the unitary S matrix S{b, s) as follows: 

'^out — S^T^inS = Tiin + T^classicali (6) 

where 

T^ciassicai = j d'^x'A{x - x'; iJ,)j{s, 6; x') . (7) 

The S matrix following from such a classical source is still an operator in the 
space of pions. Inclusion of isospin requires S{s, b) to be also a matrix in the isospace 
of the leading particles. 

The coherent production of isovector clusters of pions is described by the fol- 
lowing S matrix: 

S{s,b) = J (fe\ e)D{J;s,b){e\, (8) 

where | e ) represents the isospin-state vector of the two-leading-particle system. 
The quantity D(J; s, b) is the unitary coherent-state displacement operator defined 
as 

D{J]s,b) ^exp[J2 / dqJc{s,b]q)aJ{q) - H.c], (9) 

where aj{q) is the creation operator of a cluster c and the summation is over 
all clusters. The clusters are assumed to decay independently into c— 1,2,... pions 
and outside the region of strong interactions. 

Clusters decaying into two or more pions simulate a short-range correlation 

between pions. They need not be wcll-dcfincd resonances. The more pions in a 
cluster, the larger the correlation effect expected. 

If the conservation of isospin is a global property of the colliding system, then 
Jc{s,b;q) is of the form 

— * — * — * 

Jc(s, b; q) = Jc(s, 6; g)e, (10) 
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where e is a fixed unit vector in isospace independent of q. The global conservation 
of isospin thus introduces the long-range correlation between the emitted pions. 
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3 Distribution of Pions in Isospace 



If the isospin of two incoming particles is //a , then the initial-state vector of the 
pion field is S{s, b) \ II3), where | II3) is a vacuum state with no pions. The n-pion 
production amplitude is 

iTn{s, b;qi... Qn) = '^sil'I'^; Qi-.-Qnl S{s, b) I 1/3), (11) 

where I'l'^ denotes isostate of the outgoing leading particles. The unnormahzed 
probability distribution of producing n+yr"*", n_7r-, and no7r° pions is defined as 

W{n+nno, I'I's, Ih) = / d%dq^dq2 ...dq^l {I'hn+nn^ \ S(s, b) \ Ih) |2j(12) 
where n — n+ + n_ + no- 

Assuming further that all (/', Jg) are produced with equal probability, we can 
sum over all possible isospin states of the outgoing leading particles to obtain 



^ , , Y.ri'W{n+n_noJ'IkJh) 

Pihin+njno) = — 1^ — 13 

En+n_no E/'/^ W{n+njnQ, rr^] Ih) 



as our basic relation for calculating various pion distributions, pion multiplicities, 
and pion correlations between definite charge combinations. 

In order to obtain some more detailed results for multiplicity distributions and 
correlations, one should have an explicit form for the source function Jc(s, b;q),c = 
1,2,... . 

We shall analyze the isospin structure of our model in the so called grey-disk 
model in which 

n,{s,b)^n,{s)e{bo{s)-b), (14) 

where nc(s) denotes the mean number of clusters of the type c; bo{s) is related to 
the total inelastic cross section, and ^ is a step function. 
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4 Correlation between Neutral Pions 



We assume that pions are produced both directly and through isovector clusters of 
the p-type. 

In order to calculate the first two moments of the joint probability distribution 
-f/Zsl^O)^-) ) "we define the generating function Giig{z,n_) 

Giis{z,n_)= ^ P773(n+,n_,no)^"°, (15) 

no,n+ 



from which we calculate 



{no)n_^ ^lnGn3{'^,n_), (16) 



fln_--^JriGn,{l,n_). (17) 



The form of the generating function Gii^{z,n_)is the following 

0,,M = (/ + />x I Pl'M P (18) 

where 

2A{z, x) = (1 - x^)n^ + z{l - x^)np + 2a;% (19) 

and 

2B{z, x) = n^(l + - 2zx^) + 71^(2 - z{l - x^)). (20) 

Here denotes the average number of directly produced pions, and rip denotes 
the average number of p-type clusters which decay into two short-range correlated 
pions. The function Pj^{x) denotes the associate Legendre polinomial. Note that 

A{l,x) = B{l,x). 

The total number of emitted pions is 

(n) = + 2np. (21) 
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In Fig. 1 we show the behaviour of (no)„_ for different combination of ijiT^,np) 
when / = /s = 1. For Centauro-type behavior to appear the slope of (no)n_ should 
be negative. However this is only possible if rTp = 0. Recent estimate of the ratio of 
p- mesons to pions at accelerator energies is rip — O.lOri^. 

In Fig. 2 we show the behaviour of the dispersion D{nQ)n_ which is related to 

f2,n_ = D{no)l_ - {no)n_ (22) 

again for different pairs of {n.jr,np) and I — I3 — 1. It should be interesting to 
measure /2,„_ as it is a sensitive quatity of the pairing properties of the pions. 
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5 Conclusion 



The results of the present analysis have shown that the emission of isovector clusters 
of pions in the framework of an unitary eikonal model with global conservation of 
isospin supresses the large isospin fluctuations of the Centauro-type for pions. This 
might suggest that Centauro-type effect, if it exists, will probably appear only in 
very limited regions of phase space where isovector clusters should be missing. How 
it is possible dynamically is not clear to the authors [13-16]. 
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Figure captions : 

Fig. 1. The average number of neutral pions as a function of the number of 
negative pions for I — I3 — 1. and + rTp = 18. The curves represent different 
combinations of {nT,,np), the average number of directly produced pions and the 
average number of p-type clusters, respectively. 

Fig. 2. The correlation function for two neutral pions as a function of the 
number of negative pions for I = = 1. The curves represent different combinations 
of (n„,np). 
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